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Abstract 

We consider a coexistence of two axisymmetric liquid bridges LBi and LBm of two immiscible 
liquids i and m which are immersed in a third liquid (or gas) e and trapped between two smooth 
solid bodies with axisymmetric surfaces Si, S 2 and free contact lines. Evolution of liquid bridges 
allows two different conhgurations of LBi and LBm with multiple (five or three) interfaces of non¬ 
smooth shape. We formulate a variational problem with volume constraints and present its governing 
equations supplemented by boundary conditions. We hnd a universal relationship between curvature 
of the interfaces and discuss the Young relation at the singular curve where all liquids meet together. 
Keywords: Isoperimetric problem, multiple liquid bridges, the vectorial Young relation 
PACS 2006: Primary - 76B45, Secondary - 53A10 


1 Introduction 

Consider an evolution of two liquid bridges LBi and LBm of immiscible liquids, i (inner) and m (inter¬ 
mediate), trapped between two axisymmetric smooth solid bodies with surfaces 5 i, S2 in such a way that 
LBi is immersed into LBm and the latter is immersed into the e (external) liquid (or gas) which occupies 
the rest of the space between the two bodies (see Figure [T^). When liquid m begins to evaporate then 
LBm reduces in volume (and width). Depending on the relationships between the contact angles of both 
liquids on Si and S2 there are two scenarios for connectivity breakage of the liquid bridge m between 
the two solids. The first scenario (five interfaces) occurs when LBm splits into two parts each supported 
by a different solid (see Figure [TJ3). The second scenario (three interfaces) occurs when LBm is left as a 
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Figure 1 : a) The meridional section of two Und interfaces before LB^ rupture, b) Five Und interfaces of 
different curvatures for three immiscible liquids trapped between two smooth solid bodies with free BC. 
The endpoints Ci, (72, Cs, C4 have one degree of freedom: the upper and lower endpoints are running 
along Si and S2, respectively. The endpoints C5, Cq have two degrees of freedom and are located on two 
singular curves Li, L2, respectively, which are passing transversely to the plane of Figure. 


whole but has support only on the upper (see Figure [2J3) or lower solid. 




(b) 


Figure 2 : a) Two Und interfaces before LBm rupture, b) Three Und interfaces of different curvatures 
for three immiscible liquids trapped between two smooth solid bodies with free BCs. The endpoints 
Ci,C2, Ci have one degree of freedom while C3 has two degrees and is located on a singular curve L 
which is passing transversely to the plane of Figure. . 


Both scenarios lead to a new phenomenon which has not been discussed in literature before, namely, 
an existence of multiple LBs with non-smooth interfaces. In contrast to the known LBs with fixed and free 
contact line (CL), here one of CLs appears as a line where three interfaces with different curvatures meet 
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together. From a mathematical standpoint this singular curve is governed by transversality conditions 
(in physics they are referred to as the Young relations), and coincidence conditions, i.e., three interfaces 
always intersect at one single curve. We derive a relationship combining the constant mean curvatures of 
three different interfaces and give the interfaces consistency rules for their coexistence. Another important 
result is the vectorial Young relation at the triple point which is located on a singular curve. 

2 Variational problem for five interfaces 

Consider a functional E[r, z] of surface energy 


E[r,z] = ^ / Ejd(l)j 

j=i 

ri^i . r^\ 

+ + 

POO 

hi 



ri’2 . /■V’i 

POO 



/ + / A™dV’2 + 

/ , a^.V’2, 

(2.1) 


Jo Jij)^ 

hi 


VI 

VI 

1—1 

+ 

A\^=l\^Ra^R'i + ZS, a = 

= 1 , 2 , l = i,m,e 

(2.2) 


where r' = drj/d(j)j, Zj = dzj/d(j)j, R'^ = dRa/dil^a and Z'^ = dZa/dij^a- Throughout the paper 
the Latin and Greek indices enumerate the interfaces and solid surfaces, respectively. The surface tension 
coefficients 71 = 75, 72 = 74 and 73 denote tension at the e - m , m-i and e-i liquid interfaces, respectively, 
while 7 *^ stand for surface tension coefficients at the solid-liquid, Sa-l, interfaces (see Figure [T})). 

Two other functionals Vi [r, z] and Vm [r, z] for volumes of liquids i and m read 

/■?^i f ^2 r^l r<t>l r'l’l 

Vm[r,z\= / Vi#i - / V 2 # 2 + / y5d4>5- / V 4#4 - / Bs^dljji + / Bs 2 (iV' 2 , 

J(j>\ J(j>\ Jcpl J(j>l Jil}l dV’t 

r<p2 r^l r<t>\ r^2 

Vi[r,z\= / V2d())2 + / V3(i(/)3 + / V4#4 - / Bs^dipiV / Bs2dV’2, ( 2 . 3 ) 

where 

V, = l<j< 5 , Bs^ = ^Z'^Rl a = 1 , 2 . 

The isoperimetric problem requires to find a sef of functions fj{(l)j), Zj{(l)j), providing a local mini¬ 
mum of E[r, z] wifh fwo consfrainfs Vi[r, z] = Vi and Vm[T, z] = Vm imposed on fhe volumes of liquids 
i and m. Consider a composife funcfional 

iy[r, z] = E[r, z] - XiVm[r, z] - X^Vir, z], ( 2 . 4 ) 
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with two Lagrange multipliers Xj and represent it in the following form, 


W[r,z] = ^ / 

3 = 1 -^^. 




rvi rwi r°o 

+ / GldV'i + / + / Gld^i - 

Jo Ji^i J'4’i 

ri’i ri ’2 r°° 

/ Gldi>2 - / G^diP2 - / Gi#2, 
Jo Jip^ Jipl 


(2.5) 


where Vj = F(rj, Zj, r'-, z') and G^(i?Q, Za, R'^, Z'^) are given as follows, 


Ei-AiVi, 

F 2 — E 2 — A 2 V 2 , F 3 

— E 3 — A 3 V 3 , 

LU 

II 

LL 

-^ 3^51 + A*^, 

Gr = AiB,,+A-, 

A 2 = 

A 3 — Ai, 

A 3 B.,-a:^, 

G^ = AiB,,-A-, 

re _ \e 

'^2 - 

Gi = A^ 


(2.6) 


Si ■ 


Calculate first variation of W when the functions rj{(j)j) and Zj{(l)j) are perturbed by Uj{(j)j) and 
Vj{4>j), respectively, 

d(^j + {G\ - Gf) 6 ^Pl + {GT - Gl)6^Pl-{Gl - G^) .5V^|-(G^ - Gl)6^Pl (2.7) 




G1 -GT = A* - A™ + A 2 B,, G- - G- = A™ - A^^ + AiB,,, 


g*2-g^ = a:^-a™ + A2B,„ 

^ ^ dVj dVj , dv. , , 


51 ‘51 

e _ \m \e 


G^-Gl = A--A^^i + AiB,„ 

dVi 


( 2 . 8 ) 


The functions Uj{(l)j) and Vj{(j)j) may be derived using a requirement that the upper free endpoints of the 
first and second interfaces at Figure [Dr are running along Si and the lower free endpoints of the forth and 
fifth interfaces - along S 2 , 


rj{(tf-) = Raii’i), + Uj{(j)^) = Railpi + (JV'a), = 


' 3 y^3 


M3 


■3\^j 


dRo 

d'lpa 

dZa 




Zj{(Pj) = ZaiK), Zji<Pj) + = Zc.i'l/Ji + 6^i), Vj{(p^) = 

where a = 1 for j = 1, 2 and a = 2 for j = 4, 5. Substitute (12.8b into (12.7b and integrate by parts 


in' = E 

i=i 


, 1 



(5F,- 


<5F, 


Srj Szj 




dfi 




“i— + + 1 


dz' 


+ 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


(G1 - Gf) <5V^? + (Gr - G?) - (G*2 - G^) 5 ^ 1:1 - (G^ - G|) 

where 6 ¥/6yj = dVjdyj — d/dx /dy'^, denotes the variational derivative for the functional 
/ F ^x,yj,y'j^ dx. Since Uj{4>) and fj(<?i) are independent functions, vanishing of the integral part of 
5W in (12.9b gives rise to the Young-Laplace equations (YLE) |T|, 






Setting the remaining terms (12.111) to zero gives rise to the four transversality relations, 


^ (’^i) fr ('#'1) + ^ M) ^ M) + cr M) - G5 (I’D = 0 

Wi ^ ^ 

^ ('/’D fr (^>1) + ^ (<*8 (0l) + G 5 (VI) - G? (vl) = 0 , 

w, M + Si S + = 5 * (Vi) - Gi (Vi) = 0 . 


( 2 . 13 ) 


and one more transversality relation 


5Fi , 5Fi , . 2 ^ 5F2 , ,, 2 ^ 9^2 r.u ^Fs . ^F., 


(^?) (^ 5 ) (V)) i (-3 (Vi) 


(2.14) 


-3 (Vi) ||v .3 (Vi) (VI) (VI) I?-.. (Vi) |i -.3 (Vi) § = 0, 


In case of one liquid bridge LBm and two immiscible liquids m and e between two smooth solids 5i, S 2 
the first and forth relations in ( 12 . 131 ) coincide with those derived in |T], formula ( 2 . 15 ), while the rest 


of relations disappear. Regarding condition (12.141) . the perturbations 


and v-i 


are related in 


such a way that the three disturbed interfaces 1, 2,3 (and other three 3,4, 5) always intersect at one point. 


Ui = U 2 = U3 {(pl) , Vi {(pfj = V 2 {(pl) = V3 {(pi) , (2.15) 

(</'!) = Ui {(Pl) = U5 {(Pl) , V3 {^Pl) = V4 {(pl) = V5 {(pl) . 


Combine (12.141) . (12.151) and use independence of ui {(pf), t'l {(pi)-: U 3 {(P3), V 3 {(pf) and obtain four rela¬ 
tions. 


(^i) - II - 0 , II (,?) + II (,i) - II ^ 0 , 

|| (^i) -(<^i) - S (<^ 4 ) -(<^^) = 0 . 


( 2 . 16 ) 


dr'. 


Or! 


94 


54 


54 


Boundary conditions (BC) (12.13112431) have to be supplemented by condition of coincidence of interfaces 
in C 5 , Cq located on singular curves Li, L 2 , respectively. 


n {(pl) = r2 {(pj) = ra {(pl) , r4 {(pl) = rs ((/ig) = rs {(pf) , (2.17) 

Zi {(pl) = Z 2 {'P 2 ) = {(Pl) , Z4 {(p\) = Z^ {(p^^ = Z 3 {(pf) , 

while the angular coordinates cp^ and ipit are related by 

zi{(p\)=Zi{pj\) , ri {(p\) = Ri {(pl) , Z2{(pl)=Zi{'tpf) , r2 {(pl) = Ri (ipl) , ( 2 . 18 ) 

Z 4 (<^ 1 ) = ^2 (V’^) > ^4 {(pl) = R2 (V'a) 3 ^5 (<^5) = ^2 (V’2) 3 ’’5 (<^5) = R2 (- 02 ) • 
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Thus, we have 24 BC for the ten YLE (12.12b of the second order. Let us arrange them as follows, 


where 
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22 {(pl) = 23 (0^), 

r 2 {(pl)=Ri (V’f) , 
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( 2 . 21 ) 


?’4 {(Pl) = rs {p)l) , 

24 ((/)4) = Z2 [pjf) , 
aFs 


dR 2 

dip2 


^5 (</'5) = 
^5 (</'!) = 


ar' 

^ ^3 (</'!) , 

^2 i^l) , 


Ti {(PX) = n {(pl) , 

r 4 (0I) =i ?2 (V’l) > 

(^i) if {^D + GT (V’i) = 0 , 

^5 = ^3 {(pl) , 

^5 ( 01 ) =-^2 (V’i) > 


( 2 . 22 ) 


(2.23) 


cr (v^i) = cr (V’l) - G? (V’l), cr (v^?) = gi xpd - gt {^pd , 

GT (^2) = g*2 (V’2) - G^ , gt (v^i) = Gr (V’i) - Gi (v^i) ■ 


2.1 Curvature law and interface consistency rules 

Analysis of YLE (12.12b yields an important conclusion about the curvatures Hj of five interfaces. Con¬ 
sider (12.12b and recall that according to [Jj, Aj = 2'yjHj. Combining this scaling with (12.6b we arrive at 
relationships between the curvatures of three interfaces, 

'jiHi + 72^/2 = 73 H 3 , Hi = H 5 , H 2 = Hi. (2.24) 

Simple verification of (12.24b can be done in special cases. Indeed, if the liquids i and m are indistin¬ 
guishable, i.e., 71 = 73 and 72 = 0, then Hi = H 3 . On the other hand, if the liquids m and e are 
indistinguishable, i.e., 72 = 73 and 71 = 0, then H 2 = H 3 . In the case 71 = 72 = 73 7 ^ 0, we arrive at 
relation known in theory of double bubble 01 when three spherical soap surfaces meet at a contact line. 
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We can formulate strong statements on consistency of five interfaces based on relations (I2.24I) . 
Recall d] that there exists only one type, Nod~, of interfaces with negative H, while the other in¬ 
terfaces have positive H: nodoid Nod+, cylinder Cyl, unduloid Und, sphere Sph, or zero curvature, 
catenoid Cat. Denote by Mns"^ = {Nod+, Cyl, Und, Sph} a set of interfaces with H > 0 and by 
Mns^ = {Mns''', Cat, Nod~} a set of all admissible interfaces. The rules of interfaces consistency with 
curvatures Hi, H 2 , H 3 are given in Table, e.g., if the first and second interfaces are Cat and Nod^ then 
the third interface has to be also Nod~, but if the first and second interfaces are Und and Nod^ then the 
third interface may be any of six interfaces. 


Interfaces 

Mns+ 

Cat 

Nod 

Mns+ 

Mns+ 

Mns+ 

Mns=^ 

Cat 

Mns+ 

Cat 

Nod^ 

Nod^ 

Mns^ 

Nod^ 

Nod^ 


2.2 Standard parameterization and symmetric setup 


Consider non-zero curvature interfaces rj^cpj), 1 < j < 5, between two solid bodies, 

Za{' 4 ’a)}, 0 = 1 , 2 , and choose interfaces parameterization in such a way that the lower (j)^ and the upper 
(j)^ coordinates of endpoints Ci, (^ 2 , C 3 , 6*4 are located on the solid surfaces Si, S 2 and governed by BC 
while the other two points C 5 , Cq denote the triple points located on singular curves Li, L 2 where three 
different interfaces meet together. 

Following fil write the parametric expressions for the shape of such interfaces and rj{(f>j). 


2\H. 

M{(j),B) = {l + B)E 

1 -|- Bj cos (pj 


j\ 

(^,m\+{l-B)F 


1 + Bf + 2Bj cos 4> 




(2.25) 


, ^ BjSiTl(Pj 

^ 2\HAr 




2\Hj\rj 


(I'”)’ 

1 -|- Bj cos ( 


Bj sin (pj 


= 


AB 


(1 + ^) 


2 ’ 


r- + z- = l. 


(2.26) 


For all interfaces we have to find 24 unknowns: 15-1=14 interfaces parameters Hj, Bj, dj (due to (12.241) ) 
and 10 endpoint values (p]^4’j- These unknowns have to satisfy 24 BCs in (12.19H2.2^ . 

When both surfaces and S 2 are similar and the picture in Figure [T}) is symmetric w.r.t. midline 
between 5i and S 2 , then such setup reduces the problem above to six YLE (12.19112.211 ) for the first, second 
and third interfaces with twelve unknowns: 


</>i, di, d 2 , Bi, B 2 , B 3 , 


Hi, H 2 , 
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and c/)^ = TT, 2ds = —M{'k, \ and H'j, = ( 71 -ffi + ^ 2 ^ 1 )hz- This number coincides with twelve 

BCs which comprise ten BCs in (I2.20I2.21I) and two first BCs in (I2.22I) . Calculate the partial derivatives 
d^jjdr'p d^jjdz'j and write these twelve BCs, 


7iri {r'lB'i + ZiZ[) + ( 7 ™ - 

(</>?) = ^3 (^ 3 ) > 

72?’2 {r 2 R'i + 4 ^ 1 ) + (Tsi - 
r 2 (</>!) = n (^3) , 2:2 { 4 > l ) 


= ^3 (<^ 3 )) {(t>i)=Zi (V’i), n { 4 >i)=Ri (V'l), 

7 ™) Ri^R? + Z[^ + X 2 Z[^^^ = 0 , h = 

= Z 3 {(pl) , Z 2 {(pl)=Zi {'Ipf) , r 2 {4>l)=Ri {i/jj) , 


7iri4 + 72r2r2 - -/zr^r'^ = 0, (/>i = (j^l, (/)2 = (p^, (pz = (ph 

liriz[ + 72^24 - 73?’34 = ^ + -^24 - Xzrf) ■ 


'pJi, 

(2.27) 

V’l, 

(2.28) 


(2.29) 


After simplification we obtain 

71 cos e\ + 7™ - 7 I 1 = 0, 72 cos 6»? + 7 *^ - 7 ™ = 0, (2.30) 

ri{(pl) =r2[(pf) =rz[(pl) , zi {(p\) = Zi {p)\) , n {(p\) = Ri ['ipl) , 
zi ((pf) = Z2 ( 4 ) = ^3 (^3) > ^2 ( 4 ) =^i (^1) , 1^2 {(Ph) =Ri (V’f) , 

7 i 4 {(Pi) + 724 {(pf) - 734 {(Pz) = 0> 7 i 4 {(Pi) + 724 {(pf) - 734 (4) = 0> (2-31) 

where cos = {j'jR'i + 4 '^i) /\/+ Z'^ determines the contact angle 9{ of the j-th interface and 
Si- Two equalities in (12.301) give the Young relations at the points (7i, C2 on while two equalities in 
(12.311) represent the vectorial Young relation at the triple point C 5 located on a singular curve. Indeed, 
the latter equalities are the r and 2; projections of the vectorial equality for capillary forces f j at Cz in 
outward directions w.r.t. and tangential to meridional section of menisci. 


fi(C 5 ) + £ 2 (^ 5 ) + £ 3 (^ 5 ) = 0 , £ 4 ^ 5 ) = 7i { 4 (^ 5 ), z'(C5)} . (2.32) 

Finish this section with one more observation related the surface tensions 'yj and contact angles of three 
interfaces on solid surface. Bearing in mind that 73 cos Of + 7 *^ — 7 ^^ = 0, combine the last equality with 
two others in (12.301) and obtain, 

7 i cos 0\ + 72 cos of = 73 cos of. (2.33) 


2.3 Solving the BC equations (liquid bridges between two parallel plates) 

Making use of standard parametrization ( 12 . 251 ) we present below twelve BCs ( 12 . 3012 . 311 ) for twelve un¬ 
knowns (p\,(p\,(p\,(pf,(pf,di,d2,Bi,B2,B^,Hi,H2,imL way convenient for numerical calculations, 
^iBismcpf 72i?2sin(/)| _ 73i?3sin(/)j^ yi^icos^ 72^2 cos 4 _ 73-63 cos 4 

\IR\ \ih\ “ m ’ \m\ \^\ “ \ih\ ’ 
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(2.34) 


+ 2 Bi cos (j)l + B'l + 2 B 2 cos ( 1)2 + i ?2 y^l + 2 B^ cos cl)\ + S 3 


I Si 


IS2I 


IS, 


M(</. 2 ,Si) , ^ M{ct>lB2) , ^ M(</-1,S3) , ^ ^ M(7r,S3) 

+ “1 — - r.1 TT I - \- (I 2 — - ^1 J.J. I -T U3, (23 — —- 


2|Si| 2|S2 

M{(!)], Bj)-M(4,1 Bj 


2 \Hi 


= Zi 


2IS3I 
MjcPlBs) 
2IS3I 


2 S, 


— ds, 




1 + 2Bj cos 0] + S| 


2Si 




cos 4 >] + sin (j)^, tan 0 ^ 


1 -1 


j = l,2, 


where S 3 = S 171/73 + S272/73- 

The numerical optimization of the solution was done by a standard gradient descent algorithm. The 
cost function for the optimization problem was chosen to be the weighted sum of absolute values of the 
differences between the right and the left hand sides of the six first equations in (12.341) . 




(b) 


Figure 3: a) Five Und interfaces for three immiscible media {i - water, m - octane, e - air) trapped 
between two similar solid plates {Zi — Z 2 = 2) with free BC. b) Enlarged view of the vicinity of the triple 
point C 5 on the singular curve where three phases coexist. The angles between the adjacent interfaces 
read: 4 >i 2 = 36.76'^, 4>23 = 159.43°, 4 > 3 i = 163.81°. 

In Figure[3]we present the shapes of five interfaces of different curvatures for three immiscible media: 
i - water, m - octane (CgFfis, a component of petrol), e - air, trapped between two similar glass plates 
with free BCs and capillary parameters taken from iSj. The interfaces have the following parameters. 
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Interfaces 

7j, mN/m 

dj 

B3 

H3 

dj 

^7 

^3 

l(e — m) 

21.8 

19° 

0.959 

0.095 

-11.424 

184.39° 

208.74° 

2 {m — i) 

50.8 

39° 

0.778 

0.272 

-4.664 

186.405° 

216.34° 

3(e — i) 

72.8 

34.4° 

0.855 

0.218 

-5.599 

1 — ‘ 
OD 

0 

0 

188.14° 


The volumes of liquids confined inside interfaces read Vm = 4.009, Vi = 2.674. 


3 Variational problem for three interfaces 


Consider a functional E[r, z] of surface energy 


3 ^ /■</>) r'Pl . ri’l ri’l r° 

E[r,z\ = ^^ / ^jd4>j + / Al^d'ipi + / + / Al^Tp2, (3.1) 

i—lJidj JO Jipl d-i/’i JO J'tpl 


-JU^YJ II -'Si 

-/o 


and two functionals Vi [r, z] and Vm [r, z\ of volumes of the i and m liquids 
Vm[r,z]= I Vi#i - / V2#2 - / 


/■'/'i r<p 2 f'P] 

/ Vi#i - / V 2#2 - / 

J <t>\ ^2 


(3.2) 


2 /■'/'a r^i r'Pl 

Vi[r,z\ = I V2#2 + J ^ Vsdh- I / Bs2dV'2, 


7,/,^ .70 

where all integrands Ej, Vj and are defined in (12.2112.41) . 

Consider the composed functional W[r, z] = E[r, z] — XiVm[r, z] — z] and represent it in the 

following form, 


3=1 


/*^i /*V^i r '^2 

fjdcj)i+ Gld7Ai+/ GTdij3i+ GldiPi- G*2dV’2-/ _Gi#2 , (3.3) 

Jo J '4’i Jo 




'n 


where the integrands are given in (12.61) . 

Applying a similar technique as in section |2] we arrive at the first variation. 


<,11 


d^3 / dfj V' 


s'l' = T, 

i=i 

(G1 - GT) diPl + {GT - Gt) diPl - {Gi - Gi) 5iPl 


+ 


(3.4) 


This case does not allow a symmetric version and therefore is less reducible compared to the 5 in¬ 
terface case w.r.t. the number of unknowns and BC equations. This number equal fifteen: nine interface 
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parameters Hj, Bj, dj, and six endpoint values </>}, They satisfy fifteen BC equations 


6 fi _ 6 fi 
6 ri dzi 


dr 2 6 z 2 


' ^ W) + fr W) if W) = 0, 

' = ^3 ((/'a) > i^i) = ^3 (4) > 

zi {(j)l)=Zi (V^J), n ((/)}) =i?i (V>|), 




^ (^f) if + 

r2 (<^i) = ra , 

2:2 {cl)l)=Zi ('i/’i) , 


gf (V'f) H (</>^)+Gr (v^f) = 0, 

5^2 {4>D = Z3 (^a) , 

r2 {(pl)=Ri (V’f) , 


5ra fca 


ft (*11 + ft (^i) - ft (4) = 0 . 

t (■;>?) + ® ('#'1) - ft (<'i) =». 


dFi 
dz^ 

dR 2 dFs 
dip2 V ^2 ) drl 


I dZ 2 

' di >2 


3 t aFg 
2 / dzl 


, 2 \ ^ Qze 


Z 5 {( 1 * 1 ) = Z2 {( 1 * 1 ) , rs (0|) = i?2 (V’i) • 


= 0 , 


that gives 


(3.5) 


7icos6»J+7™-7|^ = 0 , 72 cos 0 ? + 7*^ - 7™ = 0 , 73 cos 6l| + 7*^ - 7!^ = 0 , (3.6) 

(</>?) + 72 r^ {(t*fj - 73^3 ((/> 3 ) = 0 , 71 zi (</)?) + 72 zi {(t*fj - 734 (4) = 0 > 

n (</>l) = ^^2 {(t>2j = n (4) > {(t*l) = Z2 {(t*l) = Z3 (4) > 

ri {(j)l)=Ri (V>f) , rs ((/> 2 ) =i?i (V’f) , (^1) = R 2 {(1*1) , (3.7) 

Zl {(j)\)=Zi {(l*\) , Z2 {(1*1) =Zi {((*1) , Z3 {(pl) = Z 2 {(pl) . 


Three equalities in (13.61) cannot be reduced to a single equality similar to (12.331) if the upper and lower 
solid bodies have different capillary properties, namely, 7 ^^ — 7 ^^ 7 ^ 7^2 “ i-C-> 

7i cos 0 \ + 72 cos 6 \ 7^ 73 cos 02- 


3.1 Solving the BC equations (liquid bridges between two parallel plates) 


Using a standard parametrization (12.251) and relation (12.241) for H 3 we present below fourteen BCs 
(13.613.71) for fourteen unknowns: </>}, 02 , 03 , , 02 ) </'!) ^ 2 , 7?i, i? 2 , T^a, Hi,H 2 in a way conve¬ 

nient for numerical calculations. 


1 111 Tf(0( i?,) / 

Bj = [cos 0} + sin 0] tan e{] , —-h dj = Zi j , j = 1, 2, 

B3 = [cos0| -hsin0|tan6l^]“0 ^^2\H3\^^ + da = ^2 (0^) , 

I j_r I .^1 -h 2i?i cos 0} -|- BI 41 + 2272 cos 02 + i?! 

-MtU!)-’ -MtW)-’ 
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7i-Bisin0^ 

1^11 

\/l + 2Bi cos (pl + 


72-62 sin 02 _ 73-63 sin (/jg 7ii?icos0f ^ 72-B2 cos 0| _ 73-B3 cos 
\Ih\ “ \Ih\ ’ \Ih\ + \Ih\ “ ^^ 3 ^ ’ 


1^11 


1 ^ 2 ^ , 

7 _ Y^l + 2 B 2 cos (p 2 + B 2 

w 


\Hi 

^ yi + 2 B 3 


m _ 

cos 03 + 


M{4,B2) , , M(0l,i?3) , , 


1^31 


M(02,i?i) , ^ _ Mi^lB 
2 |Fi| + ' 2 \H 2 \ 

of three interfaces of different 


2 |Fi, 

where = (-ffi 7 i + -ff 272 )/ 73 - In Figured we present the shapes 


C, 


C, 




(a) 


(b) 


Figure 4: a) Two Und and one Nod interfaces for three immiscible media (i - water, m - hexane, e - 
air) trapped between two (not similar) solid plates {Zi — Z 2 = 1) with free BC. b) Enlarged view of the 
vicinity of the triple point C 5 on singular curve L where three phases coexist. The angles between the 
adjacent interfaces reach the following values: <hi 2 = 46.04°, <^23 = 173.28°, 4 > 3 i = 140.69°. 

curvatures for three immiscible media: i - water, m - hexane a component of petrol), e - air, 

trapped between two plates composed of different materials (glass and glass coated with polymer film) 
with free BCs and capillary parameters taken from [51. The interfaces have the following parameters. 


Interfaces 

7j mN/m 



Hj 

dj 



l(e — m) 

18.4 

19° 

1.091 

0.229 

-2.521 

199.51° 

228.89° 

2 {m — i) 

51.1 

0 

0 

0.775 

0.379 

-3.257 

217.19° 

228.50° 

3(e — i) 

72.8 

49 ° 

0.841 

0.324 

-3.574 

169.79° 

211 . 11 ° 


The volumes of liquids confined inside inferfaces read Vm = 0.4377, Vi = 0.8940. 
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4 Conclusion 


We formulate a variational problem for coexistence of axisymmetric interfaces of three immiscible liq¬ 
uids: two of them, i and m, immersed in a third liquid (or gas) e and trapped between two smooth solid 
bodies with axisymmetric surfaces 81,82 and free contact lines. Assuming the volume constraints of 
two liquids i and m, we find the governing (Young-Laplace) equations (12.121) supplemented by boundary 
conditions and Young relation (12.131 ) on 5i, 82 and transversality relations (12.161 ) on singular curve where 
all liquids meet together. 

We consider two different cases when the problem allows the coexistence of five (section |2l) or fhree 
(section |3]) inferfaces. In fhe firsf case fhe problem is reduced solving 16 boundary conditions, 4 Young 
relations and 4 fransversalify relafions (12.19112.23]) . i.e., 24 equafions for 24 variables. In fhe second case 
fhis number is reduced subsfanfially, namely, 15 equations wifh 15 variables (13.51) including 10 boundary 
conditions, 3 Young relations and 2 fransversalify relafions. 

We derive fhe relafionship (12.241) combining fhe consfanf mean curvafures of fhree differenf inferfaces, 
e — m, m — i, e — i, and give consisfency rules for interface coexisfence fsecfion l^TI) . 

Anofher resulf is fhe vecforial Young relafion (12.321) af fhe friple poinf which is located on a singular 
curve. If has a clear physical inferprefalion as fhe balance equation of capillary forces. More imporfanfly, 
if gives a new insighf on an old assertion abouf fhe usual Young relafions (12.3013.61) af a solid/liquid/gas 
inferface refered by R. Finn |21 to T. Young: the contact angle at a solid/liquid/gas interface is a physical 
constant depending only on the materials, and in no other way on the particular conditions of problem, 
and a well known contradiction with uncompensated normal force reaction of solid (see @ and references 
therein). Indeed, being applied to the contact line of three continuous media, liquid-gas-solid, vectorial 
relation (12.321) assumes a singular deformation of solid surface if its elastic modules take finite values. 
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